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I. INTRODUCTION 



Lattice Monte Carlo simulations operate by dividing a finite volume of space and time into a grid, such that in a 
given direction, fi, we have points distance apart. The desire to obtain results free from uncontrolled finite 
volume contamination dictates that the product a^L^ be chosen to be suitably large in spatial directions (3 fm is 
often quoted for QCD). Controlling discretization effects similarly requires that a M be suitably small, but this must be 
balanced with the computational overhead that increases with L^. Reducing the dependence of simulation results on 
the lattice spacing is clearly advantageous, and so-called (Symanzik) improved actions may achieve this, permitting 
the use of coarser lattices without increasing discretization effects Q . 

In many cases, lattice results (such as hadron masses or decay constants) are obtained from the decay of correlation 
functions, C(r), over a range of temporal separations, r. It is a feature of such correlation functions that the signal to 
noise ratio decreases with increasing r, and beyond some r c measurements are dominated by statistical fluctuations. 
The precise value of r c depends upon many factors, including the operators correlated and the number of Monte Carlo 
measurements made, but appears to be relatively insensitive to the temporal lattice spacing, a t . As measurements 
can only be made for r an integer multiple of at , if the temporal lattice spacing is large compared to t c it will be hard 
to obtain an accurate picture of C(r < t c ). Improving the action does not help in this respect, and in addition, the 
inclusion of improvement in the temporal direction leads to the introduction of spurious (sometimes called 'ghost') 
poles in the gluonic propagator jlj . By not temporally improving the action we avoid this, but at the cost of increased 
discretization errors for given at- Controlling these, and the desire for increased temporal resolution of correlation 
functions, argues for the use of a small a t . 

We are thus motivated to choose a temporal lattice spacing that is smaller than the spatial, a s , and such 'anisotropic' 
lattices can be created by tuning action couplings in the temporal direction differently to those in the spatial ||. 
Anisotropic actions have already been successfully applied in many situations, including the glueball spectrum ||, the 
spectrum of excitations of the inter-quark potential p| , heavy hybrids |(| , and the fine structure of the quarkonium 
spectrum [Q. More recently, anisotropic lattices have been shown to be very successful in nonrelativistic QCD 
NRQCD) studies of two- and three-point correlators and finite momentum hadrons and semileptonic B decays 



A more widespread use has been hampered by the fact that the bare anisotropy (or aspect ratio) in the simulated 
action, x = Q-s/a-t, is not, due to quantum mechanical effects, the same as the measured value, \R- Typically, we 
wish to estimate the continuum limit ratio of a mass, M, to a given scale, K, using the lattice measurements, M (in 
units of at) and K (often derived from the static quark potential, and in units of a s ) respectively. Up to finite lattice 
spacing corrections, 

M = jt. (1) 
K X rK [ ' 
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We thus require xr> an d with a sufficiently small error that this does not represent the dominant uncertainty in the 
final estimate. \R mav be measured in Monte Carlo simulations, e.g. |Io[| , but it is an expensive calculation which 
must be repeated for every choice of the bare couplings. 

More generally, (lattice) perturbation theory may be used to calculate the renormalisation of quantities, and it is 
well known that with 'tadpole improvement' such calculations converge quickly to the measured data at simulated 
values of the gauge coupling, g 2 , Q. In this paper we obtain Z(g 2 ,x) = Xr/x to one loop for a wide range of 
anisotropics for SU(3) gauge theories in four dimensions. We focus on the Wilson action and a commonly used 
Symanzik improved action. Our results apply both to actions with and without tadpole improvement; in the Wilson 
case we cover both plaquette and (Landau) mean link improvement, whilst in the Symanzik improved case we discuss 
here only the mean link improved case. 

These is no reason, a priori, why the lattice spacing should be identical in all spatial directions, and indeed there are 
situations where we might choose this not to be the case. A typical example is where increased momentum resolution 
is desired for correlation functions of operators at finite momentum. Rather than increase the computational overhead 



by a global rescaling, a can be made smaller in one spatial direction 11 . Whilst we do not consider calculations 
for this case explicitly in this paper, when we describe the Feynman rules in Section [n], we allow for arbitrary 
anisotropic lattice structure as well as general actions. We use twisted boundary conditions as a gauge invariant 



infrared regulator. In Section |Tj we describe our calculation of the anisotropy from the dispersion relation of the 
on shell gluon propagator. We compare these one loop results to measurements of xr. from simulations, and show 
that the one loop result is accurate to within 3-4% over the range of couplings covered by Monte Carlo simulations. 



Finally, in Section IV we provide a summary of our findings and some conclusions. 



II. THE PERTURBATION THEORY 

It is useful to consider the derivation of an anisotropic lattice action from the isotropic continuum theory in two 
stages. We first obtain an action for an anisotropic continuum, which is then discretized. 

A. continuum anisotropy 

The starting point, and the fixed point of the lattice action, is a D-dimensional continuum field theory in an 
Euclidean space-time that is invariant under Lorentz transformations and hence isotropic. We may choose to change 
our measurement units in the continuum theory, and by different factors in different directions, which leads to the 
introduction of an anisotropy factor (or factors), \, into the action, being the ratio of the length units in different 
directions. Nonetheless, Ward identities (derived by considering anisotropics differing infinitesimally from unity) can 
be enforced to ensure that the underlying theory maintains the correct Lorentz invariance under renormalisation. 

We distinguish quantities in the isotropic theory from those in the anisotropic by the use of 'hats' in the former 
case. Although the original metric is 

Q a p = diag(+l,+l,...,+l), (2) 

we find it convenient to introduce the notation of covariant and contravariant indices. The contraction of a momentum 
and position |22]| , p • x = p a x a , must be invariant under rescaling, i.e. 

p a x a = P ^. (3) 

We can relate the rescaled fields to the original by the use of a set of vierbeins 

x^ — e^ a x , p^ e^ p a . (4) 

The metric in these variables is 

Q^u = e ll a e va , cj = det0 p „. (5) 

In the most general case there will be D — 1 anisotropy factors, but for the rescaling in the temporal direction only, 
as we consider in this paper, 
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and 

0^ = diag(i,l,...,l) , = -^- (7) 

\X J X 

Using these conventions, the natural position vector is covariant under rescaling, x = (x a ) = (t,x), so the momentum 
must be contravariant p = (p a ) — {E,p) such that Eqn. (||) is satisfied. The volume element is given by 

dx 



dx = d X Y[e^ = dx^=— . (8) 

M =i x 

The dimensionfull (colour) vector potential and derivatives behave as p, and so the Yang-Mills action becomes 

Sym = \J dScF al3 F a ^ = i J dx^fl F af} F af} = ±[ dx^fl ^ a $ VT F^F aT . (9) 
For the specific example above, 

Sym = \J^ {x 2 F 0l F 0l + F l3 F l3 ) . (10) 
B. the propagator 

To construct the Feynman diagrams for any action the gluon propagator must be computed. This is done for 
a given momentum by inverting the two-point gluon vertex, which for the isotropic (continuum) case is given by 
V( 2 ) a P(\t) = fc e fc e <5 Q/3 — k a hP . Before this can be done the gauge must be fixed and we add to the action a gauge fixing 
term and source, which in momentum space appear as 



dk 1 

(2^)4 2a' 



Sg.f. = / TK^l —M*) k>1 VM-*) - >(k)4(-k) . (11) 



The parameter a is the usual gauge-fixing parameter and, for example, a = 1 corresponds to Feynman gauge. In 
moving to the anisotropic theory, dk — dk/^/g, which affects functional derivatives with respect to the anisotropic 
source, J M . We can rescale to absorb this metric factor, which multiplies terms quadratic in by y/g. The two 
point function, V^ tJ/l/ (k), that we shall shortly derive from the action, will already contain this factor and the inverse 
propagator becomes 

r (<*W k ) = y(2W k ) + /g. r 12 ) 

a 

By illustration, the inverse propagator in the continuum for a = 1 has the form 

r <a=l )|W(k) = 5 ^^ kpk P = S J1 L 2 P + J2 kU (13) 

where the latter expression uses the isotropic momenta, which for a lattice theory we shall equate to the 'physical' 
ones. 

The propagator is 



To fix to Landau gauge we must be more careful. Consider the case where we wish to change the gauge from a to 7 
after inversion. Then 

r ( 7 )^ = p(a);u/ + f± _ £j W^fi . ( 15 ) 

We write 

G$ = Gi a J + ZV , (16) 
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and then satisfies 



The solution for D^ u is 



(r^ + - ^ V&) ( G l a J + d »p 



5% . (17) 



_ (a - 7 ) fc^fc, 
C. discretization 

The anisotropically formulated theory may be discretized in the normal way, and in these anisotropic units we set 
the lattice spacing, a, to be the same in each direction. 

On a cubical lattice Ad in dimension D (with x e Ad, fi = 1, . . . ,D) the gauge field is denoted U, 

U = {I7„(x) : tf„(x) e SU(N)} , (19) 

where C^(x) is associated with the link (x, n). 
We define the perturbative gauge field A by 

A = {A M (x + |e„) : A M (x + § e/1 ) G alg(SU(N))} , (20) 

with e p the lattice basis vectors, all of unit length and changing covariantly with rescaling. Expanding in the colour 
index, a, 

A ll = A^T a , [T a ,T b ] = -f abc T c , Tr(T a T b ) = -i(U, (21) 

where T a are the (anti-hcrmitian) generators of SU(N) with structure constants f a bc- It is expedient to associate the 
gauge potential with the centre of the link, and then 

tV(x)=expa/A M (x+^e„)), (22) 

where g is the bare coupling constant, and we have absorbed a factor of a into each component of A^. 
For a lattice with sites in the /i direction the momentum vector k is 

k =—(£-,... ,-^), 0<fc„<L„, fc^GZ, (23) 
and the sum over k stands for the sum over the components k^. In the limit that L M — > oo we have 

,24) 

The Fourier transform to momentum space is 

4(k) = £ e - ik -( x+ ^)^(x+±e M ), 

X 

A„(x + |e„) = ^ E ^^^(k) , (25) 

k 

where V = J\ is the number of lattice points. It is useful to re-express the position of a gauge potential as 
= 2(x + ^e M ), which is a D-dimensional vector with integer (covariant) components. 
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D. vertex functions 



To permit us to compute perturbation theory for a range of actions, we have developed an algorithmic method 
for expanding a general gauge theory action on a lattice in an appropriate form for perturbative calculations to be 
carried out. The approach follows closely the method and notation of Liischer and Weisz Q but is extended to 
accommodate, inter alia, anisotropy, fermionic actions, actions for non-relativistic heavy quarks (NRQCD) and more 
complicated definitions of the action in the purely gluonic sector. The algorithm is implemented in the Python 
programming language. For the work presented in this paper we briefly review the notation relevant to the present 
calculation and refer to p2| for further information. 

The lattice action for the pure glue sector can be written as a sum over contours 

S[U] = J2 C « ReTr ( p a(U;x)) , (26) 

which is defined in terms of the coupling constants, c a , and the P a , which are closed Wilson loops. 

The perturbative action is the expansion of S[U] as a polynomial in A and the coefficients of the monomials will 
determine the vertices of the theory. We denote this action as S[A] and following jl2| we write 



^^EtT E •■■ E A-\(k 1 )..AX(K)V r (k 1 ,^,a 1 -,...;k r , f i r ,a r ). (27) 



By a choice of units we set the lattice spacing to a ~ 1. The value of the lattice spacing in physical units is 
determined by a calculation of a physical dimensionfull quality and depends on g and hence on the renormalized 
coupling constant through the standard (3 function. Other quantities, such as the bare anisotropy, are determined by 
the coupling coefficients c a . 

The Euclidean Feynman rule for the r-point gluon vertex function is (— g r V r ), where the vertex V r can be expressed 
as @ 

K(ki,^i,ai; . . . ;k r ,^i r ,a r ) = ^ a ' CV(ai, . . . ,a r ) a ■ F r (ki,^i; . . . ;k r ,/i r ) , (28) 

a£S r 

where we symmetrize over S r , the permutation group of r objects. 

The C r are the Clebsch-Gordan coefficients which, owing to the reality of the action, are defined by 

C r (a u ...,a r )=Tr (T ai . . . T ar ) + (-l) r Tr (T a „ . . . T ai ) . (29) 

Under Z r , the subgroup of cyclic permutations and inversion, the C r have simple properties, 

a ■ C r — Xr(c) C r a € Z r , 

t v [ 1 for a a cyclic permutation. 

XrW) = if -nr t + u • • ' ( 3 °) 

I (—1) for a the inversion, v ' 

so it is useful to split the symmetrization operation into two steps 

K(ki,/xi,ai; k r ,/i r ,a r ) = a ■ C r (ai, a r ) a ■ Y" r (ki,/xi; kr,/j, r ) , 

tTeS r /Z r 

Y r = Y r , K a = E <?-Y r a - (31) 

a aeZ r 

The symmetrization over Z r is carried out within the Python vertex generation code, whereas any remaining per- 
mutations (for r > 4) must be carried out during the loop integration. is the contribution from the Wilson loop 
P a given by a sum of n r terms with the same momentum and Lorentz structure 

Y r a (k 1 , f ii;..r,K,Vr) = c a jr^ expT | (k x ■»? + ... + k r -<)") . (32) 

n=l r ' ^ ' 

The factor of (r!) normalises the symmetrization, and the dependence on the Lorentz indices fj,i,...,fi r has been 
suppressed. The prefactor of \ normalises Eqn. (|29|). The expansion of S[U] can thus be represented as a set of 
'entities' [/„; u", . . . , v™], n — 1, . . . , n r , where f n is an amplitude which is an integer for simple actions. The Python 
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code produces data files where these are appropriately labelled so that, given the Lorentz indices (/Zi, . . . ,fi r ) and 
the incoming momenta k 1; . . . , k,., the corresponding value of the r-point vertex function V r can be computed. The 
relevant Feynman diagrams can be constructed and the integrals over loop momenta performed either by direct 
summation over modes or using numerical integration routines. The gluon r-point functions are generated with the 
anisotropy fixed at the chosen (bare) value, and thus encoded in the amplitudes, /. We find this allows greater 
simplification of the data files produced by the Python and more efficient loop integral evaluation code. The time 
taken to rerun the vertex generation code for different \ is negligible, especially when offset against this. A more 
complete description of the implementation may be found in |13| , |l 4f| . 

The gluonic propagator is derived as per the continuum theory, using the two point vertex for the particular action, 
V^k, [i] — k, i>), and pairs of forward and backward nearest neighbour difference operators, 

A Q = (Ao.A,-) A /1 = (^,A i ), A" = (xA 0> A i ), 

X 

A 2 = A±A^ = At AJ + J2 Af . (33) 

i 

to replace the position space derivatives in the gauge fixing term. The difference operators are 

A+ = /(x + A)-/(x) 

A" = /(x)-/(x-A) (34) 

for some /(x). The net effect is merely to replace momentum components, p^, prior to any raising of the index, by 
p^ = 2 sin(p M /2) in Eqns. 



E. Faddeev— Popov ghosts 



The Faddeev-Popov ghost term is of the form S g h ost — J <ix fj a M ab ri b , where fj and rj are the usual adjoint anti- 
commuting ghost fields. The ghost fields are not observable, and form only internal lines in Feynman diagrams. We 
are thus free to choose the normalisation of the fields such that explicit factors of ^/g do not appear in the momentum 
space Feynman rules ]23| for the ghosts. The anisotropy then only appears implicitly in the raising of indices. 

The Faddeev-Popov matrix A 06 is determined by the gauge fixing condition corresponding to the choice of gauge 
in the propagator. The gauge-fixing is done by introducing the identity in the form 



1 = J da JUdA^ 8(F[A fl ] - A) exp(— A 2 /2a) 



(35) 



We use the linear gauge function (i^A^] = (A Y A^ and, as is well known, the matrix A ab is independent of a in 
this case. We denote the gauge transformation field by 



g(x) = exp(w(x)) , 
where w(x) = w Q (x)T a . For infinitesimal to the gauge field transforms as 



(36) 



ALU a 



ft M = -A+w 



n+ = (2 + a; 



J2 9 n a n (A;r b -n 



n=2 



(37) 



where A^ is the adjoint representation for the gauge field and a n is the coefficient of x n in the expansion of 



b(x) 



n=0 



cx n x 



(38) 



The Faddeev-Popov matrix is then 



8F\A l f l a ] 
M [A^] = 1 r! 



5oj b 
(A-f 



UJ = 



-A+5 ab - 



9 -Af{2 + A+) 



OO 

E 

n^2 



n _ / /\n\ab a + 
9 Un(A^) A^ 



(39) 
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The inverse ghost propagator is given by the O(g ) term and is 

A ab (x, y )=-<5 ab A 2 (x,y), (40) 

giving rise to the standard 1 /p 2 ghost propagator. 

The one-gluon vertex is given by the O(g) term which gives the contribution to the action 

^ r fcc (A -)^(l + A +/2)^) ; (41) 

and, integrating by parts, we get 

- g((A+ffj a )f abc Alv c ~ ± 5 ((A+) V)f &C (A^ C ) • (42) 

The first term generates the standard three-point vertex that one expects from the continuum but the second term is 
a lattice artifact which is suppressed by a power of the lattice spacing a as we should expect. 

At order g 2 there is a two-gluon vertex which, using o?2 = 1/12, can be read from Eqn. ( |39"| ) to be 

£((A+)V) K A l f daC f ECb ( A > b ) ■ (43) 
Higher order vertices follow a similar pattern but do not contribute to the one loop calculation we are considering. 



F. the Haar measure 



The field measure in the function integral is the Haar measure for integration over the lattice fields J7 M (x) which 
take values in the Lie group. For the perturbative calculation this measure is re-expressed as the measure for the fields 
A^x), which take values in the Lie algebra, times a Jacobian which can be expanded perturbatively and included 
as counter terms in the perturbative action. We relate the infinitesimal vectors d£ = d{; a T a and dA = dA a T a in the 
fundamental representation of the Lie algebra by 

e gd£ - — e~ gA e 9(A+tL4) (44) 

from which we derive the relation 

dA a = b(gA) ab d£ b (45) 

where again A is in the adjoint representation and the function b(x) is defined in Eqn. ( |38| ) above. The Haar measure 
is 

The Jacobian then leads to the term in the action 

/oo „ 
dx £ Tr [log 6(5^ (x))] = £ g" f3 n / dx]T Tr (l«(x)) . (47) 

Noting that 

-x-^-logb(x) = b(x) , (48) 
ax 

we find that j3 n = a n /n with a n the coefficients in the expansion of b(x) defined in Eqn. (pq) . 
The 0(g 2 ) vertex from the measure relevant to the one loop calculation is then 



\J dxJ2 ■ ( 49 ) 
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G. twisted boundary conditions 



We follow Liischer and Weisz and use twisted periodic boundary conditions for the gauge field. There is then no zero 
mode and hence no concomitant infrared divergences in the gluon self energy whilst gauge invariance is maintained. 
We briefly review these boundary conditions and refer to jL2] for further details. 

For an orthogonal twist the twisted boundary condition for gauge fields is 

C/ M (x + L v e u ) = fi^x)^ 1 , (50) 
where the twist matrices f2„ are constant SU(N) matrices which satisfy 

and z^v = exp(2irin^/N) is an element of the centre of SU(N) with = 0, 1, . . . , N — 1. The particular boundary 
conditions imposed are uniquely specified by the antisymmetric integer tensor and a complete discussion may be 
found in Jl5| . The condition that the twist be orthogonal is that n ulJ n ul j = modN, where — e fJi ^ ap n l7 p. The 
gauge potential A p (x) also satisfies the periodicity condition in Eqn. (|50|). 

Following |p^| we choose n pv — everywhere, save n\2 = —n%\ = 1. f2i and U,2 are then determined up to a unitary 
SU(N) transformation and = H4 = 1. In the case of orthogonal twist O3, can be expressed in terms of fii and 
0,2 once the values of n M! , are given. This will affect the momentum spectrum in the 3, 4-directions, but the Feynman 
rules given below are unchanged. 

The lattice is here taken to be continuous in the 3, 4-directions and of extent L sites in the 1, 2-directions. The 
momentum spectrum, k, is then continuous in k%, k^ and discrete in kx, with 

k = — - + — -, »=1,2, 0<fe i <L i , 0<m<N, k u n t eZ, (52) 
1 V Li 

with ni — n 2 — excluded to eliminate the zero mode and impose a gauge-invariant infrared cutoff momentum of 
2it/NL. Negative momentum in these directions is — fcj = (— ki,~ni), adding appropriate multiples of L and N to 
remain in the ranges defined above. 

The Fourier expansion of a colour field <j>(x) is 



(53) 



where 0n(k) is a scalar field and the sum over ki, k2 signifies the sum is over k\, k2 and the twist vector n = (ni, 122). 
The N 2 — 1 SU(N) matrices, T(n), are given in terms of fix, Q2 by 

r(n) = 2(«i+«2)(ni+"2-l)/2 fi -n2 ni ) (54) 

where z = exp(2iir/N) is an element of the centre of SU(N). We do not need to construct the T(n) explicitly but 
only evaluate the trace algebra associated with the perturbative vertices. We introduce the coefficients 

l na = Tr(rt (n)T a ) , (55) 

for which we have the relations 

Jjlnalam ~ ^nm j ( Z.R\ 

T a = ±T{n) lna , [bb) 
r(n) - 2T a7 L . 

In addition, for an adjoint field we can define the set of scalar fields labelled by n 

0™ (x) = 7na<Pa => 0(x) = —T(n)(p n (x) . (57) 

Using Eqn. ( |53"| ) ra (x) has Fourier transform 0„(k). Note that the related scalar field 

0„(x) = e - 2m( - niXl+n2X2 ^ NL (j) n (x) , (58) 
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is periodic on the lattice and has a momentum spectrum defined by the fcj in Eqn. ( |52| ) which allows the numerical 
Fourier transform to be easily computed. 

Defining the symmetric and antisymmetric products of twist vectors 

(n, m) = Uimi + n 2 m 2 + (ni + mi)(n 2 + to 2 ) , 

(n,m) = n\mi~n2m\ , (59) 
the Clebsch-Gordan coefficients given in Eqn. (E9h are modified to become 



Cj w (n 1 ,...,n r ) = ^(TT(T(n 1 )...T(n r )) + (-l) r Tr(T(n r ) . . . r(m))) . (60) 

The can be evaluated using the relations 

T(n) = 1 n = mod TV , 

Tr(r(n)) = n/0 mod N , 

T(n)t = «-i(».»)r(-n) , ^ 

r(n')T(n) = zl««'.«>-("'^))r(n' + n) , 

where (n mod N) is understood to apply to each component, rti^, and the argument of T is evaluated mod TV. We 
then derive the useful result 

C? w (-n u -n r ) = (-1)^-4 EU(«i.*N) C f w *(n lj . . . ,n r ) , (62) 

For the inverse propagator we have 

C™(n,-n) = 2z^ n V 2 . (63) 
The r-point vertex function is then given in a similar form to that in Eqn. (^fj) by 

K(ki,/ii; k r ,/i r ) = cr • Cj w (ni, . . . ,n r ) a ■ Y r (ki,fif, . . . ; k r ,/i r ) . (64) 

aeSr/z,. 

Note that the momentum arguments k implicitly define the associated twist integers n. To simplify the notation, 
we replace in most subsequent expressions the twist vector with its 'parent' momentum, understanding that only the 
twist vector will contribute in functions such as T(k). 

The structure of the vertex functions Y r is unaffected by the choice of boundary condition which is manifested only 
in the momentum spectrum used. 

A simplifying feature is to note that all diagrams contributing to an r-point function carry the same overall phase 
factor from the centre of the gauge group. These phases can be taken out as overall factors and the remaining parts 
of the Clebsch-Gordan coefficients and propagator are real. The overall phase can be restored at the end of the 
calculation. 

For the 3-point vertex on the left hand side of Fig. |l](a) we have (using momentum conservation p + k + q = 0) 0] 
Cj w (123) = lTr(r(-q-p)[r(q),r(p)] 

= ^Tr(r(q + p)r(-q - p)) (*3«<i.p>-(<J>p)) - a 4«P.q>-(P.«i))^ 

= 2isin(y(q lP )) Z T((q.q)+(p.p)+(q+p.q+p)). (65) 

In the diagram contributing to the gluon self energy of Fig. [jj(a) the phase factors for both vertices are identical and, 
with the phases of internal lines, yields an overall phase of ^sCP'P) as expected for a term in the self energy. 

For the 4-point vertex loop of Fig. |f|(b) there are three contributions from the permutations S4/Z4 corresponding 
to (1234), (1243), (1324) and in each case the Y± are real. We then find the Clebsch-Gordan factors to be 

CJ W (12M) = C 4 TH/ (1243) = 2z( k > k )/ 2 z(P>P)/ 2 

(66) 

CJ W (132A) = 2cos (f (k,p))z^ k )/ 2 z(P'P)/ 2 , 

The z( k ' k )/ 2 cancels the z~( k > k )/ 2 from the internal propagator giving an overall factor of z^- p ' p ^ 2 which is the expected 
phase. 
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H. ghost and measure Feynman rules 

We choose the Fourier representation for the anti-ghost field fj given in Eqn. ([53]) but use the conjugate twist 
matrices r^(— n). In this case the ghost propagator is real and given by 

A(p,-p) = ^. (67) 



Now consider the relevant part of the r\Ar\ vertex, 

J dxf abc fj a A b r] c = J dx2Tr(fj[A,r}]) 



2 



Tr(rt(-k)[r(p),r(q)]) J dxffcA^. (68) 



where terms such as 77k refer to the decomposition in Eqn. (pf\j, with the momentum subscript restricted to the twist 
vector. The Clebsch-Gordan structure is 

F 3 (k,p,q) =Tr(rt(-k)[r(p),r(q)]) = z-^ k - k )C^(k,p,q) (69) 

From the full structure of the f\Ar\ vertex in Eqn. ( [42| ) the momentum space Feynman vertex with momentum 
assignment shown in Fig. |l](c) is 

-i5A"co S (^)F 3 (k,p,q) , (70) 

where there is no implied sum over /i. Note that Cj (k, p, q) plays the role of a structure constant. With a real 
ghost propagator this vertex contributes the same phase factor as does the 3-point gluon vertex and so it can be 
absorbed into the overall phase factor of the diagram. 

The two-gluon vertex, Eqn. (|43|), can similarly be analysed. Assigning the ghost momentum (— k in Fig. 0(d)) as 
q, and that of the second gluon (— p in Fig. [j](d)) as r, the Clebsch-Gordan factor is 

F 4 (k, p, q, r) = J2 z ~ m) P> l ) C I W (l> ~ l ) 

I 

+C£ w (W)C? W (q,P,-Q) • (71) 
In the above we sum over a twist vector, I. The corresponding Feynman vertex is, with no sum implied over /i, 

-^V^(k-P-q-r) • (72) 

The Lorentz index of the second gluon is v. It is not so easy to extract an overall centre phase for twisted boundary 
conditions, but for the contribution to the propagator self energy in Fig. 0(d) we use Eqn. ( |62| ) to find the vertex 
contribution is proportional to z2(p>p) ICj^k, p, l)\ 2 (where we sum over twist vector I), which carries the phase 
appropriate to the self energy. 

The measure creates an insertion in the gluon propagator, and at leading order in g 2 , is 

,Z5(P.P) 

9 <W , (73) 

which carries the correct phase. As expected, these expressions are independent of the anisotropy. 

Although we do not utilise them in this paper, for completeness we also give the Feynman rules for untwisted 
boundary conditions in our notation. The ghost propagator, Eqn. (67), gains an extra factor of S a b for the ghost 
colour indices, and the measure insertion becomes 

at leading order. In the vertices of Eqns. ( |70| , |72|), the Clebsch-Gordan factors are replaced by 

F 3 (k, p, q) = 2Tr(T k [T p , T q ]) = 2C 3 (k, p, q) = -/ kpq , 
F 4 (k,p,q,r) = 4^(C 3 (k,p,e)C7 3 (q,r,e) + C 3 (k,r,e)C 3 (q,p,e)) , (75) 

e 

where, in F4, we sum over a colour index, e = 1 . . . TV 2 — 1. In both expressions, each momentum factor is understood 
to be replaced by the colour index associated with that leg of the vertex. 
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III. ANISOTROPY RENORMALISATION 



In a lattice simulation, the renormalized anisotropy is typically determined by comparing the correlation lengths 
of an operator measured along different lattice axes. In a perturbative calculation there is a much smaller choice of 
quantities sensitive to the anisotropy. 

The use of twisted boundary conditions provides one such, by providing a gauge invariant gluon mass. The 
renormalized anisotropy can be derived from the calculation of the on-shell dispersion relation for the gluon propagator 
defined in the theory with twisted boundary conditions. The details of the theory are fully discussed by Luscher and 
Weisz in jL2j and we shall follow their notation. We use a lattice of size L in the 1, 2-directions to which the twisted 
boundary conditions apply, and of extent / — > oo in the 0, 3-directions. We consider the gluon mode (called the 
A- meson in fL2[ ) which has (Euclidean) momentum 

p = (iE ,m ,0,p 3 ) , m = 2n/NL , (76) 

where Eo and p3 are continuous. In this section we understand all momentum components to be measured in the same 
units, i.e. to refer to the isotropic coor dinat e system. For clarity of presentation, however, we omit the carets used 



to distinguish such quantities in Section II A. Infrared divergences arc regulated by finite L, and mo is the pole mass 
of the bare gluon propagator. In |l2| the one loop renormalisation of mo is calculated and is used to determine the 
radiative corrections to parameters in the improved action. We follow a similar procedure and for anisotropic actions 
calculate the pole energy of the propagator as a function of p^. For sufficiently small p 3 the infrared dispersion relation 
so derived can be fitted to the standard quadratic form using the renormalized mass and renormalized anisotropy as 
parameters. 



To carry out the calculation we use Feynman gauge as described in Section II B with a = 1, but we verified that 
the results were unchanged when other gauges are chosen, as demanded by gauge invariance. The diagrams that 
contribute to the one loop gluon self energy ^^(p) are shown in Fig. [j] and the Feynman integrals were constructed 
using the vertices and rules of the previous section. 

At tree level the on-shell dispersion relation is given by r^" _1 ' ) (p) = for p of the form above, and where F^-^p) 
is defined in Eqn. jl^). In the continuum, Eqn. (fl3|), this gives 

VSPv-P" = - (~X 2 Eo 2 + P3 2 + m 2 ) = (77) 
X 

and the bare mass is defined by Eq = m>o/x at ps = 0. On the lattice for very small p% this becomes 

-(-X 2 F t (Eo)+P3 2 + F s (m )) =0 (78) 
X 

where F s ^ are determined by the details of the inverse propagator. For the Wilson action F s (m) = 4sin 2 (m/2), and 
a more complicated function with the same continuum limit for the Symanzik improved case. For the actions we 
consider below, the need to avoid extra ghost poles in the gluonic propagator (not to be confused with ghost fields) 
means that the temporal function is always unimproved: F t (E) — 4 sinh 2 (_E/2). 
Adding one loop corrections, the on-shell condition becomes 

C^r^P) = r^ =1) (^fl>™o,0,p 3 ) _ .9%, y (z£ ,mo,0,p 3 ) = . (79) 

Since we are working to 0(<? 2 ), it is sufficient to evaluate E AlJ/ (p) at the tree level on-shell energy, Eq, as per Eqn. (|7 
In general this requires taking into account the full matrix structure of Tj%) but for the form of the momentum chosen, 
Eqn. ([76]), it can be shown that of the elements Tr^J and , only is non-zero and thus that p is a zero of 

this on-diagonal element and no diagonalization of T^) is required. For given L,p3, we determine the bare pole value 
Eq by numerical solution. 

As described in |0] the field theory for the A-meson is a 2D theory, and we can write by analogy with Eqn. ( |78[ ) an 
effective dispersion relation for the infrared modes (i.e. small P3) in terms of a renormalized mass, ma, and anisotropy, 
XR- 

XR 2 F t (E R ) = p 3 2 + F s (m R ). (80) 

Using Eqn. (f79|), we have 

- (- X 2 F t (E R ) +p 3 2 + F s (m )) - g 2 Z(iE ,m ,0,P3) = 0. (81) 
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Substituting in Eqn. ( |80| ) gives 

Y 2 

- ^2 (P3 2 + F a (m fl )) + J5 3 2 + f,(mo) = g 2 xZ(iE Q , m ,0,Ps)- (82) 

We define 

Z(g 2 , x) = Xr/x = 1 + V(x)g 2 + 0(g 4 ) , (83) 
and m/j = m + g 2 mi + 0(g ), and at one loop obtain the relation 

2 V (p 3 2 + F s (m )) -mi^(mo) = x£(i£ , mo, 0,Ps). (84) 

For L -> oo we have 

X 2 / X 2 2 x 



r/(^^ _™=^.. ( 8 5) 



The values of rj(L) and m\(L) are determined by a (very good) straight line fit to Eqn. (Q). 

A. The calculation 

We work with the D = 3 + 1 SU(3) gauge theory and consider the Wilson action (W) and the Symanzik improved 
action (SI) defined in 

1. the Wilson action 
The Wilson action has a two dimensional coupling space, and is 



Sw (P, X) = -0 ( X E P M + - E ( 86 ) 

where s, s' run over spatial links in different, positive directions, P S]S ' and P Stt are plaquettes and \ is the (unrenor- 
malized) anisotropy as per Eqn. ([lO]). Spatial and temporal tadpole improvement factors, u s ,ti arising from favourite 
self-consistency conditions may be written in, but this amounts merely to a reparametrization of the same action 

_(_*__ 2"), x .*m, , 87 , 

g z J yus^ut go 2 J Ut 

and the measured anisotropy, \Ri i s invariant. We shall evaluate Eqn. ( |S3| ) perturbatively. In doing so, we span the 
full space of couplings and our calculation of Z(g 2 , \) applies equally well to any form of the Wilson action. 

To one loop the diagrams for self energy function for the gluon are shown in Fig. |l](a-e). The calculations were 
done on lattices with 8L 2 I 2 twisted momentum grid points. The loop momentum k to be summed over is 

(fc ,fci,fc2,fc 3 ) = iY no, 3L , 3L Tl2, l~ n3 ' 

0<n 0) n 3 <I, < ni ,n 2 <3L . (88) 

There is a pole in the integrands in k of the graphs in Fig. |l|(a,c) for on-shell external momentum so the kg integration 
contour is shifted by ko — > fco — iEq/2. In addition, we use the change of variables k — > k' suggested by Liischer and 
Weisz H 

K = k n ~ a v sin (M > (89) 

which gives an integrand with much broader peaks, which is easier to evaluate numerically. It is easy to see that 
a reasonable choice of parameter is a M ~ 1 — (xL^) -1 and this was found to work well, significantly reducing the 
dependence on L. 
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The integrals were done as direct summations over the discrete momentum modes but in principle they can be 
calculated using an adaptive Monte-Carlo integrator even for finite L, I . A common example of such an importance 
sampling integration package is Vegas (l6) (and described in |ljj), but this expects the integrand to be a continuous 
function of its arguments. This can be achieved by a change of variable using a stepped function Z(x, L) defined by 

Z{x, L) = int(xL) , < x < 1 , (90) 

so that Z{x, L) takes integer values in the range 0, . . . , L — 1 and a given discrete momentum component is k(x) — 
2ttZ(x, L)/L. It turns out that for one loop integrations the summation over modes is much more efficient. For two 
loops, however, it will be necessary to use Vegas to get a result of acceptable accuracy. 

We considered lattice sizes in the range 4 < L < 32 and / large enough for the error due to finite / to be essentially 
undetectable. In practice, / > 50 sufficed. We calculated the self energy at each L for a number of very small values of 
P3 and using Eqn. ( p4| ) calculated the parameters r](L) and mi(L). All computations were done on a single processor 
PC and took between 2 and 16 hours per integral, depending on L. 

A fit as a function of L allows the L — > oo values to be deduced. In Fig. |2| we show the L-dependence of rj{L) for 
X = 2, and the (excellent) fit shown is given by 

rj{L) = 0.085303 - 0.0188471-^ - 0.0073487 log ^ ^ . (91) 

For x = 1 we found that, as expected, rj(L) = for all L and that the mass renormalisation parameter mi was 

^± = -0.037923(9) , (92) 
m 

a value which agrees with that given by Snippe using the background field method |Tsj ]. Extrapolated values of r\ 
shown in Table | agree closely with results obtained by Perez and van Baal [^9|, also using background field gauge. 
For instance, we find 

r,( x = 2) = 0.085303 (1) , jj(x = 4) = 0.127901 (2) , (93) 

compared with Perez's and van Baal's values of 0.0853037 and 0.1278990 respectively. Combining our data with that 
in (nj , we find a close phenomenological fit is 

T](x) = 0.1687 (2) - 0.16397 (4)- - 0.005245 (2)4j. (94) 

X X 

2. the Symanzik improved action 
The general form of the action has, in addition to the couplings (3 and %, an additional parameter, v. 



A x,s>s' v J 

where R ss , s ' and R ss j are 2x1 loops. The anisotropy renormalisation factor, Z(x, /3,v), is a function of all three 
parameters. 

It is usual to restrict the simulated parameters to a two dimensional surface in the coupling space. The particular 
surface is determined by the chosen method of tadpole improvement, demanding v = u s {g 2 , x) f° r the reparametriza- 



tion in Eqn. (87). If, for given anisotropy, u s has a perturbative expansion u s {g ) — 1 + d s g + 0(g ), then we can 



rewrite the action as 

Ssi(P,X,v) = S SI ((3, x ,v = l) + g 2 AS SI + 0(g i ) 



ASsi = -0d. Y -l\ x Rss,t + -(Rss, s > +R s >s>.s] 



(96) 



This form is numerically more convenient than the reparametrized form, in that the dependence of the action on the 
gauge coupling is lessened. 
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In this paper we calculate Z(g 2 ,x) over two such surfaces: v — 1, corresponding to an action without tadpole 
improvement (in which case AS si makes no contribution), and the surface v — u s corresponding to Landau mean 
link improvement where u s t are defined as the expectation values of the traced link matrices in the Landau gauge 

u 

u s = (U S ) L , u t = (U t ) L . (97) 

U s j are gauge fields in the spatial and temporal directions. The actual values of u St t are established self-consistently 
by an appropriate iteration procedure, usually for fixed /?o and \o (which we comment upon in Section^). As AS si 
is 0(g 2 ) relative to other terms in the actions, we treat it as a counterterm, giving rise to a gluon propagator insertion 
shown in Fig. |l|(f). The coupling, d s , is determined numerically in each case for an identical lattice jl3], Q. 

In both cases, the actual calculation proceeds as for the Wilson action (the only changes in the loop integration code 
being the PYTHON-generated input files). In Fig. |^ we plot r)(x) as a function of x- The results are also presented in 
Table |. The slightly larger errors in the full calculation reflects a greater uncertainty in the extrapolation in L due 
to reduced data. 

These values can be fitted extremely well, as is demonstrated in Fig. |[ by 

n(x) = 0.0602 (1) - 0.0656 (2)- - 0.0237 (1W (98) 

X X 

without the mean link counterterm contribution, and 

riix) = 0.0955 (4) - 0.0702 (16)- - 0.0399 (14)-=, (99) 

X X 

when this is included. 



B. comparison with measured data 



The goal of this paper is to provide determinations of the renormalied anisotropy using perturbation theory. It is 
clear that for sufficiently small couplings the one loop result will be sufficient, but at couplings more typical of current 
lattice simulations we must check what systematic truncation errors are introduced by neglecting contributions of 
higher orders. We can do this by comparing the perturbative predictions with what measurements have been made 
in Monte Carlo simulations. 

We start with the Symanzik improved action, for which measurements of xr. have been published in |], derived 
from torelon dispersion relations and the so-called sideways potential. 

We have obtained as a perturbative series, Z(g 2 ,x, v ) = 1 + vix, 1 *)^ + 0(g 4 ). Such a series, however, shows 
poor convergence and thus large truncation errors. We can re-express it as a power series in the 'self-consistent', or 
'boosted' coupling, go 2 : 

Z(go 2 ,x) = 1 + V(X, v = u s )g 2 + O( 5o 4 ) , (100) 

using the definitions in Eqn. (|S7|). For appropriate choices of tadpole improvement schemes, the convergence is 
much improved. In resumming the series, we formally require only that u s> t are self-consistent for fixed /3ojXo u P 
to a suitable order of perturbation theory, and that v — u s to the same order. There are arguments, however, for 
satisfying these conditions at a numerical level for greater convergence. The one loop result is, of course, unchanged 
under rcsummation. 

The action in || [T(| contains explicit factors of u S}t and before applying Z it is necessary to rescale xo to x using 
the simulated values of the mean link factors. In Table || we compare the one loop perturbative determination with 
the published measurements. The couplings are larger than those used in many simulations (giving a relatively coarse 
spatial lattice), and the anisotropy factors are at the higher end of those usually employed, making this a rather 
stringent test of the perturbative series. Despite this, we find this error to be uniformly small. For xo < 4 the error 
is 2%, and for xo = 6 it is only 3%. If we attribute the difference between the measured and one loop estimate to be 
that due to truncation of the perturbative series, it is clear that the two loop contribution to Z is very small. The 
range of couplings, however, is insufficient to allow us to reliably estimate the two loop coefficient from the data. 

A comparison of the one loop perturbative and measured anisotropy has been carried out for the Wilson action by 
Klassen j2l|]. Corroboration of this analysis is hampered by a lack of published tadpole improvement factors. Whilst 
X is quoted directly, u s> t are required to obtain the self-consistent coupling, go 2 , appropriate to simulated values of 
fi. We may estimate the tadpole factors, but in doing so we introduce a second systematic error into the estimate 
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of xr- This arises from the difference between our estimate and the numerically self-consistent tadpole values, ft 
is in addition to the systematic error arising from the one loop truncation of the perturbative expansion of Z. We 
anticipate the latter being of order 3%, and we thus require values of u s , t at least this accurate. Such an estimate 
can come from a two loop perturbative estimate Oj. Using this we have confirmed using Landau mean link 
improvement Klassen's analysis (which used plaquette tadpole factors). We find the combined two loop tadpole, one 
loop anisotropy prediction for \R to be correct to within 5% for (3 > 5.6, and to 2% for (3 > 6.3 for 1.5 < % < 6. 
With the correct u Stt values, a more precise estimate of the truncation errors in Z(xt9 2 ) would be possible, and we 
expect them to show that the one loop anisotropy renormalisation calculation is as accurate for the Wilson action as 
we have demonstrated it to be for the Symanzik improved case. 

IV. SUMMARY 

We have presented perturbative calculations of the renormalisation of the anisotropy in simple gluonic lattice 
actions for SU(3) in D = 3 + 1. These have been carried out by studying the infrared dispersion relation of the 
on-shell gluonic propagator, regulated through the use of twisted boundary conditions as in fl2[| . We have reviewed 
the derivation of the Feynman rules for general actions and anisotropy structure, using both twisted and untwisted 
boundary conditions. We have (briefly) outlined our method for generating the gluon r-point vertex functions using 
a Python code and evaluating the loop integrals using a compiled programming language. The strength of this 
dual-pronged approach is the generality; at the price of not carrying out optimisation specific to a particular action, 
a flexibility to rapidly change actions is acquired. Further details of the implementation will be given in |l4[ . 

We have focussed our attention on two actions, the Wilson plaquette action and an action (Symanzik) improved 
at tree level through the addition of 1 x 2 loops extending two lattice spacings in spatial directions only. This lack of 
symmetry is necessary to avoid spurious ghost poles in the gluonic propagator, and in consequence even an isotropic 
lattice will not remain so under quantum corrections. Temporal improvement is not crucial, of course, when a t is 
small. 

Both actions are considered with and without tadpole improvement of the links. In the case of the Wilson action 
no explicit reference to u Syt need be made in the perturbation theory. Specific tadpole improvement schemes manifest 
only in rescaling the gauge coupling. 

For the Symanzik improved action, the lack of temporal improvement means all reference to ut can be reparametrized 
away, but dependence on u s remains in the form of a counterterm to the (tadpole) unimproved action. The presence 
of the counterterm necessarily specialized the perturbative calculation to one form of tadpole improvement, and we 
selected self-consistent Landau gauge mean links for this work. 

There are also advantages in using this reparametrization when tuning actions such as Ssi for Monte Carlo sim- 
ulation, when explicit tadpole factors u Syt are often present, reparametrizing to remove dependence on u t makes 
the finding of a self-consistent u s for fixed (/3, \) a one parameter problem. Similarly tuning u t is done by varying 
(/3o,XOjMt) together along the single parameter curve of constant (/?, \t u s)- This two stage approach is much faster 
than a simultaneous exploration of the two parameter space. 

For the Symanzik improved action, we presented results over a spread of bare anisotropies 1 < Xo < 8, and 
give interpolating fits, Eqns. (|9§| , |99| ), respectively with and without (Landau) mean link tadpole improvement. 
Comparing these numbers with measurements in simulations shows that the one loop perturbative determinations of 
the renormalized anisotropy are accurate to within 2% for \ < 4 across a range in couplings and lattice spacings much 
larger than would typically be used in lattice simulations. For % = 6, which is larger than that currently employed in 
large scale simulations, the deviation is 3%. 

In the case of the Wilson action we have verified existing results (calculated using the background field method) 
by Snippe Jl8[ and Perez and van Baal p9[. Using an interpolating fit, Eqn. (|94|), we have carried out a comparison 
with lattice simulation results by Klassen }21[ . We find the agreement to be approximately 3-4% for gauge couplings 
in the typically simulated range and \ < 6, but this uncertainty includes a systematic error from having to estimate 
tadpole improvement factors, and the error in the anisotropy alone is almost certainly less. 

We conclude that the renormalisation of the anisotropy, which is difficult to calculate accurately in simulations, is 
well described by lattice perturbation theory. If the anisotropy is to be used, for instance, to set the scale in lattice 
calculations as per Eqn. (Q), then a 3% systematic error in \R IS sufficiently small that is unlikely to represent the 
dominant uncertainty in the final estimate, and the one loop determinations presented in this paper are all that is 
required. 
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< e > (f) 

Haar measure one-loop mean link 



FIG. 1: Self energy graphs at 0(g 2 ) for the gluon propagator. Graph (f) arises from treating the spatial mean link factor as a 
counterterm, and has a coupling proportional to the one loop expression for this. 
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FIG. 2: Anisotropy renormalisation for the Wilson action as a function of L for bare anisotropy \ — 2. 



X 


Wilson 


Spatially Improved 


(a) 


(b) 


1 





-0.02907 (1) 


-0.01461 (30) 


1.5 




0.005654 (5) 


0.03163 (30) 


2 


0.085303 (1) 


0.021446 (4) 


0.05037 (30) 


2.5 




0.03027 (1) 


0.06097 (30) 


3 




0.035692 (8) 


0.06753 (30) 


3.5 




0.03955 (1) 


0.07218 (30) 


4 


0.127901 (2) 


0.04251 (1) 


0.07567 (30) 


5 




0.04851 (2) 


0.08000 (30) 


6 




0.04851 (2) 


0.08279 (30) 


8 




0.05172 (5) 


0.08637 (30) 



TABLE I: Anisotropy renormalisation, rj, for various x f° r Wilson and Symanzik improved actions. The latter is shown (a) 
without and (b) with the spatial mean link counterterm contribution of Fig. [j](f). 



19 



0o 


Xo 


Us 


u t 


-l 
a s 

(MeV) 


(mcas. ) 

Xr 


X 


nix) 


(pert.) 

Xr 


1.7 


4 


0.7164 


0.98295 


661 (11) 


3.56 (2) 


2.915 


0.0667 


3.601 




6 


0.7100 


0.99158 


779 (28) 


5.28 (2) 


4.296 


0.0770 


5.463 


1.8 


4 


0.7279 


0.98243 


797 (21) 


3.61 (2) 


2.964 


0.0673 


3.629 




6 


0.7216 


0.99208 


839 (9) 


5.31 (2) 


4.364 


0.0773 


5.489 


2.4 


3 


0.7868 


0.9771 


1200 (50) 


2.71 (3) 


2.416 


0.0596 


2.776 



TABLE II: Renormalized anisotropy as determined from torelon dispersion relations and the sideways potential [l(]], and 
from one loop perturbation theory. 
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FIG. 3: Anisotropy renormalisation for the spatially improved action as a function of the bare anisotropy x> showing fits 
quadratic in x~ l - 



